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We propose a method for controlUng a class of low temperature chemical reactions. Specifically, 
we show the hydrogen abstraction channel in the reaction of formaldehyde (H2CO) and the hydroxyl 
radical (OH) can be controlled through either the molecular state or an external electric field. We 
also outline experiments for investigating and demonstrating control over this important reaction. 
To this end, we report the first Stark deceleration of the H2CO molecule. We have decelerated a 
molecular beam of H2CO essentially to rest, producing cold molecule packets at a temperature of 
100 mK with a few million molecules in the packet at a density of ^10^ cm~^. 

PACS numbers: 82.20.-w, 33.80.Ps, 39.10.+j, 32.60.+i 



Recent exciting developments in ultracold matter re- 
search include the creation of ultracold molecules by 
magneto-association 1], leading to molecular Bose- 
Einstein condensation |2|. Despite the rich physics 
demonstrated in these systems, all are characterized by 
spherically symmetric interactions whose effect is in- 
cluded through one parameter, namely the s-wave scat- 
tering length. By contrast, the permanent electric dipole 
moment possessed by polar molecules permits long-range 
and anisotropic interactions and enables new methods 
for external control in an ultra-cold environment 3]. 
The electric dipole-dipole interaction (and control over 
it) gives rise to unique physics and chemistry including 
novel collision and chemical reaction dynamics. Lack 
of spherical symmetry in the interaction causes collid- 
ing molecules to be attracted or repelled depending on 
their relative orientation. Thus, an external electric field, 
which orients the molecules, will have a profound effect 
on the molecular interactions 4]. Furthermore, because 
the transition states of a chemical reaction often involve 
specific orientations of the molecular dipoles, an exter- 
nal electric field may shift the energy barrier to reaction, 
making a particular reaction pathway more or less favor- 
able. A Stark decelerator [3, 1^, 13 producing cold polar 
molecules with tunable and well-defined translational en- 
ergy is thus an ideal tool for the study of low (or negative) 
barrier chemical reactions. We note other experiments 
utilizing photo-association have succeeded in 

producing ultra-cold heteronuclear alkali dimers. 

In this letter, we present a study of the collision and 
reaction properties of OH with H2CO at low tempera- 
tures, followed by a detailed description of experiments 
required and carried out so far to probe these novel dy- 
namics. The calculations suggest for the first time that 
chemical reactions, as well as collision cross sections, can 
be altered orders of magnitude by varying either the 
molecular state, or external electric field strength. We 
specifically consider the H-abstraction channel in the re- 



action of H2CO and OH: H2CO ^ OH CHO + H2O. 
This reaction not only represents a key component in the 
combustion of hydrocarbons, but also plays an important 
role in atmospheric chemistry where it is the primary pro- 
cess res pon sible for the removal of the pollutant, H2CO 
m El in El El E 113. Near room temperature the 
rate of this reaction is weakly dependent on temperature, 
suggesting the barrier to the process (if any) is very low 
. Measurements of the thermal activation energy, Ea , 
are inconclusive 0, 0, 0, Q, 0, 0, , ranging from 
Ea/R = 750 K to -931 K, with the most recent mea- 
surement giving a value of -135 K 18]. Here R is the 
universal gas constant. The most accurate calculations 
19] predict the energy of the transition state for the ab- 
straction to lie between -700 K and +60 K relative to the 
reactants. 

Experimentally, both OH 0, 0, and H2CO (this 
work) molecules have been produced at low temperatures 
via Stark deceleration. In the case of spontaneous reac- 
tion {Ea < 0) m El El El, by magnetically trapping 

OH in the presence of a tunable bias electric field and 
"bombarding" the trap with decelerated H2CO packets 
the OH - H2CO total scattering and reaction rates can be 
mapped as a function of collision energy and applied elec- 
tric field. Thus Ea can be measured. If En > 0, as some 
measurements and theory suggest Q, 0, E^, the 
collision energy tuning afforded by the Stark decelerator 
provides a direct way to measure the energy barrier to 
reaction. Unlike thermal kinetics studies, which rely on 
fitting the Arrhenius formula to reaction rates, the Stark 
decelerator can be used to tune the collision energy above 
and below the threshold energy. 

As shown in Fig. Q^a), H2CO is a near-symmetric 
prolate top molecule, with nearly degenerate rotations 
about the b- and c- axes. The dipole moment, /i^, along 
the a- axis causes these nearly degenerate, opposite parity 
states to experience a large Stark shift with an applied 
electric field, as shown in Fig. ^b-d). Here, the states 
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FIG. 1: (a) (Color online) H2CO. (b-c) Stark shifts of the 
low- lying rotational states of H2CO. States of interest for this 
work are shown in bold. 



are labeled by their zero-field identity in the \JrTnj > 
basis, where r = Kg — Kc {Ki is the projection of J 
along the ith axis) [2l|. Of particular interest to this 
work is the |lil > state (Fig. dc)), which is the upper 
component of the lowest J level of ortho-formaldehyde 
and is hence well populated in our supersonic expansion. 
This state experiences a large Stark shift (1.32 cm~^ at 
125 kV/cm), making it excellent for Stark deceleration. 
There are other states accessible via, for example, stim- 
ulated Raman adiabatic passage that offer an improved 
Stark deceleration efficiency (|222 >) or a good candidate 
for implementing an AC Stark trap (|OoO >) ji^ . 

OH-H2CO Theory. The Hamiltonian, H, used in 
the scattering calculations takes the following form: 

H = T -\- HoH + HH2CO + H Stark + Hdd + ^sr + Hchem- 

(1) 

Here T is the kinetic energy, Hqh and Hh^co are the 
Hamiltonians of the separated molecules, and H stark de- 
scribes the action of the electric field on the individual 
molecules. Hdd represents the long-range dipole-dipole 
interaction between the molecules, which depends on the 
orientations of the molecules relative to the field axis, as 
well as the distance R between them. This will be ex- 
pressed in a basis of molecular eigenstates (including the 
effect of the field), as was done in Ref. . In addition, to 
prevent the dipole-dipole interaction from overwhelming 
the short-range interaction, we replace its varia- 
tion with 1/{R^ + Odd) for an appropriate constant Cdd 
[2^ . The chemical reaction dynamics are modeled in a 
schematic way by the terms Hsr and Hchem- The goal of 
Hsr is to mimic an anisotropic interaction with a barrier, 
not necessarily to characterize the OH-H2CO system. 
This approximation is justified for the low reaction bar- 
rier where the dipole-dipole interaction dominates over 
internal molecular dynamics. Thus, we construct a po- 
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FIG. 2: (a) Computed elastic and inelastic OH - H2CO colli- 
sion cross sections as a function of external electric field at 1 
mK. (b) Chemical reaction cross sections from the schematic 
H-abstraction model. 



tential as an expansion into Legendre functions, 
H,r = Vo{R)-2V2{R)C2o{eA). 



(2) 



Here (^, (j)) represent the spherical angles of the relative 
coordinate between the molecules, referred to the electric 
field axis, and C20 is a reduced spherical harmonic ji^ . 
The isotropic part of this potential carries the barrier at 
a radius R ^ R^ ^ 
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The height of the barrier depending on Db is fully ad- 
justable with Wb characterizing the range. In the follow- 
ing, we set the barrier height equal to the threshold en- 
ergy of the incoming molecule. The radial dependence of 
V2 is chosen to have a magnitude comparable to Vq^ but 
to change sign as R crosses R^j. Thus the polar molecules 
attract in a head-to-tail configuration for R > R^j, but 
attract in a tail-to-tail configuration upon crossing the 
barrier. Finally, chemical reactivity {Hchem) is modeled 
by coupling to a purely repulsive artificial channel. The 
coupling to this channel is represented by a decaying ex- 
ponential in R, which is of negligible size for R> R^. To 
complete a chemical reaction, the molecules must cross 
over the barrier, reverse their relative orientation, and 
find the artificial channel. 

Collision cross sections for OH-H2CO scattering in this 
model are shown in Fig. El^a). Here both molecules 
are assumed to be in spin- stretched, weak- field-seeking 
states: \FmF-,pcirity >= |22,/ > for OH, and |lil > 
for H2CO. The cross section for elastic (spin-changing) 
collisions is shown as a solid (dashed) line. The cross sec- 
tions are large and comparable, and exhibit strong mod- 
ulations as a function of an applied electric field. This 
has been explained for OH-OH scattering on the basis of 
purely long-range scattering by the dipole-dipole interac- 
tion Thus, these results are robust and independent 
of the detailed form of the short-range interaction. 

The sensitive dependence of cold molecules on inter- 
nal state and electric field engenders new prospects for 
probing and controlling the chemical reaction itself. To 



illustrate this, we show in Fig. [^fb) the chemical re- 
action cross section versus electric field. The collision 
energy is assumed to be 1 mK. (Note that, because the 
reaction mechanism is schematic in our model we report 
only relative values for the cross sections.) The solid line 
is for collisions in which both OH and H2CO molecules 
are in their strong-field-seeking states. At this energy 
the field has a profound infiuence on the reaction cross 
section, mostly by accessing a very large number of reso- 
nant states. Considering that each resonant state spans a 
different region of configuration space, mapping out res- 
onances such as these would serve as an extremely sensi- 
tive probe to transition states. The dashed trace in Fig. 
Eljb) shows the reaction rate for weak-field-seeking states 
of the reactants. In addition to probing many resonances, 
this cross section rises by nearly four orders of magnitude 
as the field increases. This rise is also a direct conse- 
quence of the dipole-dipole interaction. In Refs. 0,11^ it 
was shown that for polar molecules in weak-field-seeking 
states, the dipole-dipole interaction can "shield" these 
molecules from getting close enough together to react 
chemically. Thus the cross section is strongly suppressed 
relative to that for strong-field-seekers. However, as the 
field is increased, the inner turning point of the relevant 
potential curve moves to smaller making the shielding 
less effective, and hence, the reaction more likely. 

While the above discussion is relevant for future studies 
in a low temperature trap environment, another impor- 
tant capability for current experiments is the direct con- 
trol of the H-abstraction reaction barrier height through 
the application of an external electric field. The basic 
H-abstraction reaction mechanism, follo wing Ref. [l^ . 
is shown in Fig. |2l We note that Ref. [2^ presents a 
slightly different H-abstraction configuration, but with 
the same end-to-end dipole coupling scheme. During the 
H-abstraction process, the energy required to rotate one 
dipole moment (OH) versus the other (H2CO) enables 
our proposed control of reaction dynamics using an ex- 
ternal electric field. Hence, the two abstraction configu- 
rations present no difference other than the magnitude of 
the shift we can create in the reaction barrier. The impor- 
tant characteristic to note is that the hydrogen-bonded 
complex (HBC) in the second panel (Fig. forms along 
an attractive direction of the two electric dipoles, but 
in the transition state (TS) the OH dipole has essentially 
fiipped its orientation relative to the H2CO dipole. While 
an external electric field, which orients the molecules 
would allow (and perhaps even encourage) the forma- 
tion of the HBC; it would add an energy barrier to the 
formation of the TS, and thus the H-abstraction chan- 
nel through the OH dipole-field interaction. Based on 
the bond geometries shown in Ref. 14], the addition to 
the energy barrier in going from the HBC to the TS is 
^ 1.8/io//|^|, where j^qh is the expectation value of the 
dipole for the OH molecule. For a high, yet attainable, 
electric field of 250 kV/cm, the additional barrier energy 



Reactants 



HBC 



t 

i i ' 




Products 

t 



6 

i 



FIG. 3: (Color online) The reaction mechanism for the H- 
abstraction channel following Ref [ij- 



would be 10 K. Thus, assuming a reaction barrier at a 
H2CO colhsion speed of 223 m/s (60 K 19]) the applied 
electric field would shift the required collision velocity to 
234 m/s, well within the Stark decelerator resolution. 

Experiment. The decelerator used for H2CO is sim- 
ilar to the apparatus in our previously described OH ex- 
periments 0,0,123. Molecules in a skimmed, pulsed su- 
personic beam are focused by an electrostatic hexapole 
field to provide transverse coupling into the Stark decel- 
erator. The Stark decelerator is constructed of 143 slow- 
ing stages spaced 5.461 mm apart with each stage com- 
prised of two cylindrical electrodes of diameter 3.175 mm 
separated axially by 5.175 mm and oppositely biased at 
high voltage (± 12.5 kV). Successive stages are oriented 
at 90^ to each other to provide transverse guiding of the 
molecular beam. The geometry of the slowing stages pro- 
vides an electric field maximum between the electrodes 
with the field decreasing away from the electrode cen- 
ter. Switching the electric field when the molecules are 
directly between two adjacent stages (no net decelera- 
tion) is denoted by a synchronous phase angle of (j)o = 0^ 
(often referred to as bunching), while switching the elec- 
tric field when the molecules are between the electrodes 
(maximum deceleration) is denoted by (j)o = 90^. 

The H2CO molecules are produced from the cracking 
of the formaldehyde polymer to produce the monomer, 
which is passed through a double u-tube apparatus [2^. 
Xenon at 2 bar pressure is fiowed over the collected 
H2CO, held at 196 K where H2CO has -2.7 kPa (20 
Torr) vapor pressure. This Xe/H2C0 mixture is ex- 
panded through a solenoid type supersonic valve produc- 
ing a beam with a mean speed of 350 m/s with approxi- 
mately a 10% velocity spread. 

H2CO is detected using laser-induced fluorescence. 
The molecules are excited from the |lil > ground state 
by photons at 353 nm generated from a frequency- 
doubled, pulsed-dye laser pumped by NdiYAG laser to 
the A^A2 electronically excited state with one quan- 
tum in the out-of-plane bending vibrational mode. 
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FIG. 4: Stark deceleration of H2CO molecules, showing (a) 
bunching (^o = 0"^) and Stark deceleration at (j)o = 30° (offset 
down for clarity), (b) Stark decelerated packets of H2CO 
molecules at the output of the decelerator. 



Approximately 40% of the excited H2CO decays non- 
radiatively 29], while the remaining molecules emit dis- 
tributed fluorescence from 353 nm to 610 nm (ioj . 

Shown in Fig. ^are time-of- flight (ToF) measurements 
taken at the exit of the decelerator, which is ~0.9 m 
down- stream from the valve. In these figures the time 
axis is relative to the beginning of deceleration. In Fig. 
Wi^) where the decelerator is operated at = 0^, a large 
peak in the ToF curve denotes the arrival (marked by an 
arrow) of the phase stable packet whose density is ~10^ 
cm~^ and contains a few million molecules. It is impor- 
tant to note that even though the speed of this packet 
is 340 m/s, 10 m/s lower than the pulse mean speed, no 
deceleration has occurred. Instead the molecules in the 
pulse whose speed was in the range from 340 m/s ± 25 
m/s were loaded into the packet as detailed in our earlier 
work [23 • Also shown in Fig. ^a), but offset for clarity 
is the ToF curve for deceleration at (j)o = 30^. Here the 
arrow denotes the arrival of a slowed H2CO packet with 
a mean velocity of 204 m/s. Displayed in Fig. ^b) are 
ToF curves for the slowed molecular packets as a function 
of (f)o' The minimum speed shown is 26 m/s for = 50^, 
resulting in a temperature of ~100 mK. 

To observe the OH-H2CO collision and reaction dy- 
namics requires monitoring both the OH population and 
the production of the formyl radical (CHO). We have 
developed sensitive fluorescence imaging techniques for 
OH I3, and because CHO has a low ionization poten- 
tial (8.1 eV), high signal-to- noise ratio photo- ionization 
detection of the reaction can be implemented. Neverthe- 
less, given the current available densities of the react ants 
the observation of the dynamics proposed here will be 
daunting. However, further improvements such as multi- 
ple trap loading and secondary cooling of the molecular 
samples should make these observations possible. 

In summary, we have detailed new methods to study 
interspecies molecular collision and reaction dynamics, 
and have shown how control of these processes may be 
achieved. We have also demonstrated preparation of both 
cold polar molecules required for these studies. Exper- 
iments are currently poised to begin exploring the rich 



physics presented by the dipolar interaction. 
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